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REPORT  SUMMARY 


The  objective  of  this  project  is  to  investigate  the  effects  of  using  an 
electron  beam  to  artificially  disturb  the  upper  atmosphere,  thus  simulating 
natural  disturbances  such  as  auroras.  This  rej>ort  involves  both  a  general 
study  of  the  effects  of  electron  radiation  on  the  upper  atmosphere  and  the 
specific  application  of  these  effects  to  an  experiment  performed  bv  the 
National  Aeronautic  and  Space  Administration. 

The  pur|)ose  of  this  paper  is  to  present  an  analysis  of  electron  beam 
interaction  with  the  upper  atmosphere.  In  order  to  better  understand  the 
phenomenology’  that  would  be  involved  when  the  atmosphere  is  artificially 
excited  by  an  electron  beam,  specific  aspects  of  an  experiment  performed  by- 
Hess  et  al.  (1970)  are  analyzed.  Electron  beam  spreading,  energy-  deposition, 
and  penetration  are  studied,  Reactions  involving  air  chemistry  and  radiation 
which  result  from  this  interaction  are  also  studied.  Information  from  the 
NASA  "artificial  aurora"  experiment  has  been  analyzed  and  compared  to  the 
results  of  this  study. 

The  study-  is  reported  in  three  parts:  the  spreading  and  penetration 
of  an  electron  beam  in  the  upper  atmosphere,  together  with  the  effects  of  the 
motion  of  the  beam  accelerator  are  discussed  in  Section  3.  The  atmospheric 
chemistry  that  is  involved  due  to  the  electron  interaction  is  discussed  in 
Section  4.  Finally,  these  calculations  are  compared  with  NASA  observations 
in  Section  5.  While  the  NASA  data  was  quite  limited,  the  description  pre¬ 
sented  of  the  electron  irradiation  of  the  upper  atmosphere  is  in  good  agreement 
with  these  data. 

The  width  of  the  electron  cloud  as  well  as  the  energy  deposition  has 
been  adequately  modeled  by  this  analysis.  The  effects  of  the  geomagnetic  field 


ix 


strength  and  orientation  has  been  Included.  The  altitudes  at  which  electrons 
deposit  most  of  their  energy  Is  a  function  of  the  Initial  electron  energy.  This 
penetration  has  been  modeled  by  taking  Into  account  atmospheric  density  varia¬ 
tion  with  altitude,  tts<ng  this  formulation,  calculations  of  the  energy  deposited 
as  a  function  of  altitude  have  been  made  which  show  good  agreement  with 
observed  data 

The  chemistry  tnd  radiation  Involved  due  to  electron  Irnullatlon  has 
been  studied  In  detail.  A  complete  description  of  the  visible  emission  Involved 
is  presented.  This  calculated  emission  Is  In  good  agreement  with  that  observed, 
both  In  Intensity  and  In  temporal  behavior.  With  this  agreement,  confidence 
has  been  gained  for  predicting  the  emission  Intensities  during  both  natural  and 
mnn-mado  disturbances  of  the  upper  stmosphere. 
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5!  CTIOS  I 


INTRODUCTION 


t'Rwt  atmospheric  (»f»fiK«rwi  during  quiescent  and  itliturM 
condition*  arc  prrwrtl)-  umSrr  Intensive  »lwt)  ({Knihue,  19711, 

atrfatror  (Stair  ind  Qravtn,  1967)  and  rwidhnroe  (Saiikr  and  Narvtsi,  |970| 
(datlorma  (atr  nMtlnnd  I"  - Mttot  measurements  of  pltvttua  awl  |on*deitsllv 
profile*  (On*  hue  ct  •!.,  1179.  Iridmar  e|  al. .  It* la),  concent ra lion*  of 
r^lnor  specie*  aod  mclasiahlr*  (7 ||.(  ct  al.,  It7fa.  awl  radiation  from  Ibese 
a|*mra  (llaMha?,  IW*i.  l*i»rltn*ot»  r|  al  ,  ItToli).  Theoretical  calculation* 
ha*r  *  ho*  n  that  ion1  concentration  measurement*  dwrlnc  qulrKrnl  condition* 
arc  In  r*rcllcnt  agreement  »nh  straih  Male  chemistry  awl  radiation  calcula¬ 
tion*  (limahue,  196C,  t  l»tc\,  IW*.  titMhir,  I9.*ai.  I|«n*r%cr.  much  lc** 

I*  wndcr*food  during  almo*t  her  lc  disturbance* ,  such  a#  aurora* 

Measurement*  taken  from  ground*  and  aircraft -t»a*ed  ln*irumcnl*  o| 
ihc  radiation  emitted  during  auroral  adit  it*  hate  sfarsn  rapid  temporal  tarla 
lion*  in  l he  radian* c  intensities.  ||o*eiet,  three  measurement*  bate  leen 
unable  lo  Joc'.te  lie  precise  altitude  ol  the  radiating  specie*  (Chamberlain, 

!•*!•  INmahue,  lt<l|.  RorM dope  measurement*,  on  the  other  hand,  hate 
yielded  information  on  the  altitude  dependence  o|  the  radiating  steclr*  folder 
and  Narclsl.  l?*o  leldman  ei  al. ,  197  lat  fart  onlv  ten  limited  information  on 
on  temporal  *  aviation*  Sue  h  hate  also  hern  launched  malnlv  Into  (he  quid  |om 
break-up  stage  of  auroral  actlrttr. 

A  complimentary  technique  to  tlete  dfrrcl  measurement  program*  I* 
lo  artificially  create  an  aurond  .rpr  disturbance  a*  ha*  been  proved  fa  Sfalr 
fit” II  Using  this  Irchr.Mprr,  an  arllflrial  aurora  uould  be  produced  fa  launching 
a  high-altilude  probe  *llh  an  onfasard  r  ectron  learn  system  lo  irradiate  the  atmt**- 
pfare  »llh  electron*.  An  et)<rrlmrni  similar  lo  this  proposal  uas  performed 
by  NASA  In  It^.f  (Hr**  el  al  ,  1971.  Ilatts  el  al..  1970.  197  1».  A  rodkd 
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wmnai  mm  fJrttm  hcfl#ribi  mam  mm*  *  ir^iK  Hft  trlilt-  i*t 

W™r  inpbfbrti  i,  rt*  >  »  In  w  *j  4  «■  Win  1«»  #nafw‘  NS 
fi*M  H<*«  fw  f|ti(r>M  I>1>  «*»  dirrrlx*  l!^<  »m  ^14^1# 

ftoU  *ai  Ihf  .Irfui^  *U  <«—  nw*  a**'  Uo  *****  •* 

rttnuiiot  Iran  KM  lr  |1*  In  t«  illl'.ndr  Additional  ikuil*  rrpnilac  (fell* 
njTrimrftt  arr  ptrtralrd  ta  tbr  nril  «*rUM 

T*  |»iT*ral  aliacty  ctmtofa  a  ikorrllrtl  n»«*lr|  for  (Hr  rlrvtrnn  kim 
imrrarlK*  *tth  (Hr  nmtapHrrf  Thr  inrMr*  a  df«ffi|i»«i  of  (hr 

k'felK  l«l  catrnt  of  the  ialrnrlkM  rr{|M,  (Hr  atmo*j»hrric  Hh*Watl«» 
m*£i  iV  rkmtjiin  and  rwtiitlm  follow  lac  (Hr  rlfrltwi  irradiation  Tfe(»  mmlrl 
la  irnfwl  |a  <*r|*|j  la  5rdK«*  3  and  I  .N«*clflr  calculation*  of  (hr  aimoajdwric 
interaction  under  (few  roirtllka*  of  tHr  NA5A  rapertment  a rr  prrfnttnrd  »*t n* 
(Hi*  modrl  |*i  SrfUwi  3,  the  iatrnait)-  an*!  tlmr  ilr|r«lfnrr  of  (hr  rfitWr  rad  la 
(Ion  arc  nlniliM  and  mm|«arrd  to  (Hr  olarruKwi*  madr  in  (hr  NASA  rxpert* 
mcnta.  Particular  intcrcat  la  paid  »o  (Nr  mtchialim  for  |»n*H»clns  U|  !fe  radia 
(loo  which  *j*  found  10  (Nr  nirrlmmi  (n  |rr»l*t  aftrr  rlrciron-tirajn 
irradiation  Nad  cra«r*f 
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SECTION  2 


NASA ‘ARTIFICIAL  AURORA*  EXPERIMENT 


On  January  20,  19*9,  at  9  45  t'T  -l.’i  A  M  ,  1ST),  an  Ac  nd**e  :t'»0 
rocl.ri  carrying  an  electron  l>cnm  accelerator  «i»  launched  from  Wallop 
Island,  Virginia  The  rocket  trajectory  had  an  azimuth  97°,  reaching  an 
■i*>ltov  of  2*9  km  at  a  distance  27  km  cant  of  the  launch  |*»int.  Near  n|*ogee. 
the  rocket  was  maneuvered  no  that  the  onboard  electron  gun  directed  the  i>cnm 
downward  along  the  geomagnetic  field  The  relation  between  the  rocket  tra¬ 
jectory  and  the  magnetic  field,  which  hsd  a  local  Inclination  of  09  5°  and  a 
declination  of  -8°  I*  shown  In  Figure  I 


P* 

it  *  ■  «  »•  C  J> 


\ 


I  *- 

fji 


K  '  * 

.t 


The  electron  leant  accelerator  c*ommenced  ojh* ration  20(»  sec 
after  launch  A  21  pulse  aeifuence,  each  consisting  of  a  pulse  evert*  2.7  sec, 
w*a  rct»eated  for  252  sec.  The 
pulse  voltages  and  currents  were 
varied  from  1  to  8.7  kV  and  from 
1.5  to  490  mi,  respectively.  Ihilso 
duration  w*s  cither  0. 1  sec  or 
0.99  sec.  The  strongest  pulse 
(numbers  11  and  21  in  the  pulsing 
science)  were  operated  at  490  mi 
and  8. 7  kV  for  0.  98  sec.  Four  of 
these  pulses  were  observed  before 
the  electron  beam  accelerator 
passed  out  of  the  view  of  the  ground- 
based  detectors.  On  closer  examina¬ 
tion  of  the  data,  one  additional  trace 
operating  at  490  ma,  4.9  kV,  for 
0. 1  see  was  also  observed. 


1**4  *  *t4»*  &  HAVA 
f  lK»»Oh  m  w 

**AC*A’C0  COIUMV-J 


•AllO^S  i^kAhC  IA</HCk 

••  r*A^ai%C«Tv 


‘  r 


•  ■ 
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T#v| 
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Figure  1.  Sketch  of  the  Rleetron  Beam 
Accelerator  Trajectory  and 
the  Magnetic  Field  Orienta¬ 
tion  Described  bv  Davis 
(1971) 
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Measurements  on  the  irradiated  region  were  made  from  ground-based 
stations  located  at  ^ranklin  City  and  Igor  using  image-orthicon  systems  with  a 
framing  speed  of  30  frames  per  second.  Figure  2  shows  the  response  curve  for 
the  S-20  photocathode  image  orthicon  and  lens  system  (Davis  et  al  . ,  1970). 

The  detection  sensitivity  of  this  system  was  estimated  to  be  500  rayleighs  at 
the  |>eak  of  the  response  curve. 


| 
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Figure  2.  Relative  Response  of  the 
S-20  Photocathode  Multi¬ 
plier  and  Lens  Transmit¬ 
tance  (Taken  From  Davis 
Hal.,  1971) 


The  signals  recorded  during 
the  experiment  are  presented  in 
Table  1.  The  radiation  was 
observed  from  narrow  tubes 
extending  from-  104  to  130  km 
in  altitude.  The  measured  appar¬ 
ent  width  was  GO  to  240  meters  with 
an  average  of  about  130  meters. 

The  column  brightness  while  the 

electron  beam  was  on  was  estimated 

1 

to  be  approximately  10  to  15  kilo- 

ravleighs  after  correcting  for  instru 

mcnt  scnsitivitv  and  was  due  mainlv 
+  '  '■ 

to  N2  first  negative  radiation.  The 
total  radiation  was  observed  for 
approximately  1.3G±  0.07  sec, 

0.38  sec  longer  than  the  0.98  sec 
irradiation  time. 


TABLE  1 


OBSERVED  WIDTH  AND  UPPER  AND  LOWER  ALTITUDES  OF  SIGNATURE  PULSES 

v. 


ALTITUDE 

AT  TOP,  km 

ALTITUDE 

AT  BOTTOM,  km 

WIDTH 

meters 

FRANKLIN  CITY 

1st  Signature 

127  -  133 

105.9 

120 

2nd  Signature 

130  -  138 

105.5 

126 

3rd  Signature 

123  -  131 

113.7 

— 

4th  Signature 

137.5 

108.6 

132 

IGOR 

1st  Signature 

123.2 

104.1 

123 

2nd  Signature 

125.3 

103.9 

(64) 

3rd  Signature 

126  -  128 

103.6 

117 

4th  Signature 

128.1 

104.5 

54 
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SECTIOK  3 


THEORETICAL  MODEL  OF  ELECTRON  BEAM 
INTERACTION  WITH  THE  ATMOSPHERE 


Previous  studies  of  electron  interaction  with  the  upper  atmosphere 
have  been  performed  with  particular  emphasis  on  the  relationship  to  auroral 
activity.  Calculations  of  ionization  profiles  due  to  electron  excitation  have 
been  made  by  Chamberlain  (1961).  The  interaction  of  a  single  collimated 
electron  beam  with  the  upper  atmosphere,  which  is  analyzed  here,  presents 
a  somewhat  more  manageable  problem.  This  model,  moreover,  is  directly 
applicable  to  the  NASA  electron  beam  experiment. 

The  various  effects  of  the  electron  beam  interaction  with  the  upper 
atmosphere  can  be  modeled  by  laboratory  measurements.  These  effects  include 
atmospheric  density  variations,  the  magnetic  field,  and  chemical  and  radiative 
processes  following  t.ie  irradiation.  Experimental  data  on  electron  interaction 
with  air  has  been  obtained  in  the  laboratory  by  GrUn  (1957,  and  by  Cohn  and 
Caledonia  (1970). 

The  energy  deposited  into  the  upper  atmosphere  by  an  electron  beam 
has  been  calculated  by  Aerger  et  al.  (1970)  whose  purpose  was  to  describe 
electron  transport  from  very  high  altitudes  into  the  upper  atmosphere  as  in  the 
case  of  auroral  activity.  Using  the  Monte  Carlo  technique,  they  calculated  the 
spread  of  a  beam  of  monoenergetic  electrons,  injected  downward  from  300  km 
along  the  magnetic  field,  as  a  function  of  altitude.  A  vertical  magnetic  field 
strength  of  0.6G  was  assumed.  Under  the  conditions  of  constant  density  and 
the  absence  of  a  magnetic  field,  their  calculations  were  in  good  agreement 
with  laboratory  measurements.  A  simplified  formulation  of  this  problem  based 
on  their  results  is  presented  here  to  facilitat?  analysis  of  the  chemistry  and 
radiation. 


Preceding  page  blank 


In  the  absence  of  a  magnetic  field,  the  energy  deposition  in  air  by  an 
electron  beam  is  nearly  spherical  in  shape,  its  diameter  being  determined  bv 
both  the  electron  energy  and  air  density.  The  range  of  the  beam  measured  by 
Griin  (1957)  was  found  to  fit  the  relation 


Range  =  9.  2  x  1016  E7^4  n"1  cm 

o 

for  constant  air  density  in  the  energy  range  between  2  and  40  keV. 


(1) 

Here, 


Eq  is  the  initial  electron  energy  in  keV  and  n  is  the  number  density  per  cm°. 
For  application  to  the  atmosphere  with  varying  density,  Equation  (1)  must  be 
modified  to 


n(s)  ds  =  9. 2  x  1016 


(2) 


where  sq  is  the  position  of  the  electron  beam  accelerator,  sf  is  the  end  of  the 
electron  beam  range,  and  n(s)  is  the  atmospheric  number  density  along  the 
path  s  ,  in  our  case,  along  the  magnetic  field.  For  a  given  n(s)  ,  s^  is 
determined  numerically  using  Equation  (2). 

From  both  laboratory  measurements  (Grim,  1957;  Cohn  and  Caledonia, 
1C,  u)  and  theoretical  calculations  (Berger  et  al. ,  1970;  Spencer  and  Coyne, 
1962),  the  spatial  distribution  of  energy'  deposition  has  been  found  to  be  independ¬ 
ent  of  the  electron  beam  energy.  Tiie  energy  deposition  along  the  initial  beam 
direction,  T(z/L)  ,  is  plotted  as  a  function  of  z/L  in  Figure  3.  Here,  z  is 
the  particle  thicknesc  in  units  of  molecules  per  unit  area,  and  L  is  the  beam 
range  in  molecules  per  unit  area.  The  energy  deposited  at  negative  z/L  is 
due  to  backscattered  electrons.  The  quantity  T(z/L)  is  normalized  by  the 
relation 
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I 


/ 


T(z/L)  d(z/L)  =  1  . 


-oo 


The  normalized  particle  thickness  is  hence  defined  as 


s 

l 


n(s)  ds 


i 


n(s)  ds 


z/L  = 


/'■ 


9.2  x  1016  E^4 
o 


n(sji  ds 


(3) 


s 

o 

Thus,  for  a  constant  density  gas,  z/L  reduces  to  a  simple  dimensionless 
distance  with  L  defined  by  Equation  (1). 


Figure  3.  Energy  Deposition  as  a 

Function  of  Distance  Away 
From  the  Electron  Beam 
for  Constant  Density  and  No 
Magnetic  Field 


Consider  the  case  in  which 
the  electron  beam  is  initially  directed 
along  a  magnetic  field.  Electrons 
scattered  perpendicular  to  the  mag¬ 
netic  lines  of  force  are  turned  by  the 
field  and  the  spherical  shaped  could 
is  thus  narrowed  to  a  column.  The 
range  and  energy  distribution  along 
the  field  do  not  change  and  are  still 
represented  by  Figure  3  .  Electrons 
that  are  scatterer]  by  the  atmosphere 
receive  a  velocity  component  perpen¬ 
dicular  to  the  magnetic  field  and  thus 
make  circular  (or  spiral)  orbits 
around  the  field  lines  These  orbits 
have  a  Larmor  radius: 
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(4) 


m  cv 
e 

qB 


± 


where  q  and  are  the  electron  charge  and  mass,  respectively;  B  is  the 
magnetic  field;  c  is  the  velocity  of  light;  and  v^  is  the  component  of  velocity 
perpendicular  to  the  magnetic  field.  In  terms  of  the  electron  energy,  E  , 


R 


L 


c  sin  0 
qB 


(2me  E) 


1/2 


(5) 


where  0  is  the  angle  between  the  direction  of  the  motion  of  the  electron  and  the 
magnetic  field.  In  the  vicinity  of  Wallops  Island,  B  =  0.6  gauss,  and  RT 

Li 

becomes 


Rl  =  1.  7  sin  0  E1^2  meters  (6) 

where  E  is  expressed  in  keV.  For  the  case  in  which  the  mean  free  path 
between  collisions  is  larger  than  the  Larmor  radius,  the  lateral  spreading 
of  the  beam  is  approximated  by  the  relation 

=  2  R^  (number  of  collisions)  *^2  (7) 

where  R^  ,  the  electron  Larmor  radius  averaged  over  both  the  energy  during 
slowdown  and  the  angle  0  ,  is  approximately 

RL  =  1. 1  (Eq)  1//2  meters  .  (8) 
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The  number  of  collisions  that  an  electron  encounters  during  slowdown  is 
approximately  (Berger  et  al. ,  1971) 


(number  of  collisions) 


=  25.  6 


(9) 


Combining  Equations  (7),  (8),  and  (9),  the  lateral  spread  becomes 


S  — 
J. 


11  E^^  meters 


(10) 


Note  that  for  the  case  in  which  the  collisional  mean  free  path  is  much  larger 
than  the  Larmor  radius,  the  lateral  spread  is  independent  of  the  density.  The 
electron  cloud  when  directed  along  the  magnetic  field  becomes  a  column  with 
a  diameter  S±  at  the  widest  point.  z/L  «  0.5  . 

The  average  energy  W  required  to  form  one  ion-electron  pair 
(Landshoff  et  al. ,  1966)  is 

W  =  33.73  +  [i.53(Eo  -  E.)  1^2]  eV  per  ion  pair 

where  Eq  is  the  initial  beam  energy  in  keV,  and  E.  =  0.015  keV  is  the  mean 
ionization  energy  for  air.  The  ion-pair  formation  rate  in  the  atmosphere  by 
electron  beam  irradiation  is 


I  E  T(z/L)  , 

ion -pairs /cm  -sec  = - - -  “y  J 

IV  S  2  dS 


v-3 


'5/2. 


=  1.86  x  10  u  1Eq  '  T(z/L)  n(s) 


(12) 


11 


where  Eq  and  I  are  the  beam  voltage  in  kilovolts  and  current  in  amperes, 
respectively,  and  n(s)  is  the  atmospheric  number  density  per  cm^. 
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SECTION  4 


atmospheric  chemistry  associated 

WITH  ELECTRON  IRRADIATION 

The  chemistry  associated  with  the  electron  beam  irradiation 
determines  the  radiation  observed  from  the  ground.  Consequently,  in  order 
to  predict  or  analyze  the  observed  signal  from  electron  irradiation,  a  detailed 
study  of  the  chemistry  must  he  performed. 

The  beam  electrons  which  ionize  the  atmosphere  use  the  remainder 
of  the  35  eV  energy  yielded  in  the  creation  of  one  ion  pair  in  dissociation  and 
in  excitation  of  electronic,  vibrational,  and  rotational  modes.  This  interaction 
leads  to  reactions  which  produce  visible  radiation.  The  chemistry  that  follows 
the  electron  irradiation  may  be  described  by  coupled  second-order  differential 
equations  of  the  form 


8M 

at 


M  + 


9M 

at 


chem 


(13) 


The  first  term  on  the  right-hand  side  accounts  for  the  diffusion  effects  to  be 
discussed  at  the  end  of  this  section.  The  second  denotes  the  change  in  con¬ 
centration  of  M  due  to  chemical  reactions.  This  term  is  of  primary  interest 
in  predicting  the  electron  interaction  with  the  upper  atmosphere. 

Since  the  altitudes  at  which  the  electrons  deposit  their  energy  (for 
the  purposes  of  this  analysis)  are  above  100  km,  the  chemistry  consists  mainly 
of  binary  reactions.  Thus,  to  simplify  the  analysis,  both  radiation  and  electron 
interaction  with  the  atmosphere  are  treated  in  the  same  form.  The  general 
reaction  for  producing  species  M  is  written  as 

k. 

A.  +  B.  -i  M  +  C.  +  D.  (14) 

i  i  J  ]  ' ' 
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where  the  subscript  i  denotes  any  species  that  is  involved  in  forming  species 
M  and  k.  is  the  rate  constant  for  the  reaction.  The  reverse  reaction  is 
written  as  a  separate  equation 

k. 

Cj  +  M  -J  A.  +  B.  (15) 

with  subscript  j  denoting  a  species  involved  in  removing  M.  The  rate  constants 
k.  and  k^  are  related  by  the  equilibrium  conditions.  The  net  change  in  the  con¬ 
centrations  of  species  M  due  to  chemistry  becomes 


aM/9t|chem  =  £  ki  [A^Bj]  -  £  ^  [m][cJ  (16) 

where  the  brackets  denote  the  concentration  of  the  species  designated  within 
these  brackets.  This  binary  reaction  scheme  has  been  generalized  to  handle 
from  zero  to  three  reactants  and  from  zero  to  three  products. 

In  order  to  treat  the  electron  beam  interaction  with  the  atmosphere, 
an  artificial  chemical  reaction  with  zero  reactants  and  a  production  rate  equal 
to  that  calculated  from  Equation  (12)  was  used.  This  production  of  io.  "  was 
then  weighted  by  the  relative  production  of  each  ionic  species. 

Reactions  leading  to  ionization  and  their  relative  production  rates 

are: 
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Reaction 

Relative 

Production 

Rate 

Reference 

+ 

N„  +  e  +  e 

2  p 

0.75 

(Lnndshoff  et  ol. ,  1967) 

au> 

VN2 

-»  N+  +  N  +  e  +  e 

P 

0.  19 

(Rnppetal.,  1965) 

(H2) 

Ng  <B2  +  c  + 

0.0G 

(Borst  and  Zipf,  1970a) 

(H3> 

°2  +  e  +  ep 

0.70 

(Landshoff  et  ol. ,  19G7) 

(IM) 

V°2 

—  0  +  o+  +  e  +  e 

P 

0.29 

(Rapp  et  si.,  1965) 

atr») 

°2  (b*Ig)  +  e  +  ep 

0.01 

(Borst  and  Zipf,  1970b) 

(H0) 

and 

e  +  O 

P 

—  0+  +  e  +  e 

P 

1.0 

(Kioffer  and  Dunn,  1966) 

(H7) 

The  subscript  p  denotes  the  fast  beam  electron.  Reactions  (R3)  and  (RG) 
illustrate  examples  of  the  production  of  excited  states.  In  these  cases,  the 
ions  radiate  in  the  first  negative  bands 

N2  <®2  <>  -  N2  <X'  2g> 

and 

0>4V-°;  <x+ng) 

These  radiation  sources  are  particularly  useful  as  their  intensities  are  directly 
related  to  their  ionization  rates. 


*"8 

T  *  G.  7  x  10  sec  (RS) 

(Borst  and  Zipf,  1970a) 


r  *  1.  2  x  10  b  sec  (R9) 

(Borst  and  Zipf,  1971) 
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Tbe  electro**,  both  prtnury  and  iwkAi ry ,  »l*o  ir»*dbrc 
ditaoclatlon  and  eieitatton  ol  tbe  neutral  m*$?clc*  Tbe  important  rv*etl<m* 
include 


Re art  tor 

Relative 

FroducUo* 

JUi _ 

Reference 

• 

O,  •  r  -  O*  u  •  e 
* 

0  31 

(t-aodatioff  et  al  ,  IW7) 

(*irat  and  7.tpf.  197 1| 

fit  10} 

O  •  r  —  0(^  •  e 

0  90S 

(Henry  et  al. ,  19i9) 

(HI  1) 

O  *  e  —  0(  *  |f»  •  e 

>30 

(Henry  ct  al  ,  If*V) 

(HI2) 

N,  *  *  -  X_  *  e 

*  * 

0.A 

f*n|ar*4|  ef  al. ,  19*?} 

(HI  3} 

Nj*  r  —  X  •  X  •  e 

0,  S7 

fttlnler,  If**} 

(RH> 

The  atotnW?  auygen  produced  in  (fraction  jHIOj  i*  InyghJv  rxrifrd  rMomm*  nod 
7if»l.  1*7 1),  l»d  one  Kill  o|  the  tjxytrn  »lom*  produced  in  l(rft<rlinn  (It  10)  in 
believer!  to  be  in  the  *f*fe  (Farltnaon  et  *1, .  I»7ob,  Mumma  and  7,ipf, 

19*1.  Feldman  ct  al..  197la|  Tbe  Twin  atate*  pruriMod  in  Itearlion  (IU3> 

are  N2(C3I^),  NjII^IM,  and  SJ(A3r*» 

3  3 

Tbe  S^fC  Ij^l  and  N^cR  I^J  «tatr*  decay  rapidly  In-  radiating  the 
second  and  lira!  poalitre  band  araiem,  rcapectltrlv 

NVc3V  “  N2(,|3V  *  *>.*  »  I**"  (Rl'd 

(Andevaon.  197 1) 

In  tbe  ultraviolet,  and 

Njin'iy  -  r  -3*lo‘Sec 


(Anderaoo,  197 1 ) 


(Itir.) 


jgl 

la  the  infrared  ia d  vtaihle  »*»alnn»«ha,  fbe  N_{A  r**  *t*«c  It  long  lived  and 

I  *  « 

may  be  a  aource  of  Of  £)  Ihro#  Hm?  rwua«  f»„  ./kinaon  and  7. Ijd.  1970) 

•'‘jtA  1^)  •  O  —  N,  •  Of1!*  >  9  a  |0  **  cm^/aec  (R17) 


The  reirtlMi 

V  •  a  ■«  V  *  •  « 
'  * 

•  * 


(HIS) 


and 


O  •  r  -  Ol  •  e 
3  3 


(HI  9) 


are  t*mwfl*nl  HfenHnoil  rtrit/Hoe  reaction*  which  acctainl  (nr  I  he  accondarv 
elect  row  rnpryv  Aie»iuim»!f!»  «(  the  (»eam  energy  goe*  into  vibrational 
ear  Itaft  ion  |Chr«,  IS»M| 

The  Ion*  tirwfciced  in  HcacUon*  (HI)  through  (H7)  undergo  Ionic  reaction* 
or  recombine  »f*h  the  *{««  wrrotKtary  electrons  The  |irtnrl|«l  reaction*  ami 
their  rate  oa«*>aat*  fat  Wt  \)  are  ll»ted  )n  Tatde  I  The  relative  rale  nl 
diaacic  lathe  recomNinat ion  reaction#  and  lon'etchangr  reaction  nubatantltlly 
affect*  the  relathe  concent  rattan  of  the  ionic  #t«-clc# 
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TABLE  2 

ION  EXCHANGE  REACTIONS 


Reaction 

•Reaction  Rate 
(em3/sec) 

Reference 

X  "*  N2  +  N+ 

1  x  10‘H 

(DASA,  1967) 

(R20) 

0  -  NO "  +  N 

2.5  x  10‘10 

(Norton  and  Barth,  1970) 

(R21) 

O  —  NO+  +  N(2D) 

7  x  10‘U 

(Norton  and  Barth,  1970) 

N2  + 

o  -  n2  +  o+ 

1  x  lO-12 

(Goldan  et  al . ,  19GG) 

(H23) 

°2-N2  +  0^ 

3  x  10"10 

(Dunkin  et  al. ,  19G8) 

(R24) 

02  -  NO+  +  NO 

1  x  10‘17 

(DuiUin  et  al. ,  19G8) 

(R25) 

NO  -  N'  +  NO+ 

2 

3  x  10*10 

(Goldan  et  al.,  19GG) 

(R26) 

’  O  —  N  +  O^ 

IXIO’12 

(Schmeltokopf  et  al. ,  1968) 

(R27) 

°2  -  N  +  °2 

3  x  10"10 

(Dunkin  et  al.  ,  1969) 

(R28) 

02  -  NO+  +  O 

3  x  I0_li 

(Dut.kin  et  al. ,  1968) 

(R29) 

N+  + 

02  —  O  +  NO 

1  X  10‘:2 

(DASA,  1967) 

(R30) 

NO  —  N  +  NO+ 

8  x  10*10 

(Dunkin  et  al. ,  1968) 

(R31) 

NO  -  N2  +  O 

3  X  10"12 

(DASA,  1967) 

(R32) 

NO  -  0+  +  N2 

1  x  10"12 

(DASA,  1967) 

(R33) 

N2  -  NO+  +  N 

1  x  10‘12 

(Schmellekopf  et  al. ,  1968) 

(R34) 

o+  + 

°2  -  °  +  °2 

2  x  lO*11 

(Dunkin  et  al. ,  19F8) 

(R35) 

NO  —  O  +  NO+ 

-12 

1.3  x  10 

(Goldan  et  al. ,  1966) 

(R36) 

N  -  NO+  +  O 

1.  8  x  lO*10 

(Goldan  et  al. ,  1966) 

(R37) 

°2  + 

NO  -  NO+  +  02 

8  x  10'10 

(Goldan  et  al. ,  1966) 

(R38) 

N2  -  NO+  +  NO 

1  x  10'17 

(Feldman  et  al. ,  1971b) 

(R39) 

Rates  given  at  300 


K  unless  otherwise  noted, 
18 


TABLE  2  (Continued) 


Reaction 

Reaction  Rate 
(cmtysec) 

Reference 

N2  +  e  -  N(2D)  +  N 

2.8  x  10"7(T/300)“0,2 

(Mehr  and  Biondi,  1969) 

(R40) 

O*  +  e  -  0(1D)  +  O 

2.  1  x  10  ”7  (T/300)'0, 7 

(Feldman  et  al.,  1971a) 

(R41) 

O*  +  e  -  0(1S)  +  0 

2.2  x  L0‘8  (T/300)"0,7 

(Feldman  et  al,  197 la; 
Donahue  et  al.,  1968b) 

(R42) 

NO+  +  e  —  N(2D)  +  O 

3  x  10 "7  (T/300) ”1,0 

(DASA,  1967) 

(R43) 

NO+  +  e  N  +  O 

1  x  10_7fT/300)"1,0 

(DASA,  1967) 

(R44) 

0+  +  e  +  M— -0  +  M 

1  x  10_26( T/300)”  2,5 

(DASA,  1967) 

(R45) 

N+  +  e  +  M— -N*  +  M 

1  x  10”26(T/300)”2' 5 

1 

(DASA,  1967) 

(R46) 

) 


Radiative  and  chemical  processes  also  involve  the  neutral  species 
excited  or  formed  by  the  electron  beam: 


Reaction 

*  Reaction  Rate 
(cm*Vsec) 

Reference 

2  ^ 

N(D)+02~N0  +0 

6  x  lO-12 

(Nicholet,  1^65;  Black,  1969) 

(R47) 

N  +02(a1Ag)  -*  NO  +  O 

3  x  lO-15 

(Hunten  and  McElroy,  1968; 

Clark,  1970) 

(R48) 

n+o2~~fo  +  o 

2.4x  iu  *'•  ,69°/T 

(Kaufman,  1969) 

(k49) 

N  +  NO  -  N„  +0 

£  & 

1.5  x  10_12T1/2 

1 

(Feldman  et  al-*  1971b) 

(R50) 

0  +  NC>  -  N02  +  hy 
(green  continuum) 

6. 4  x  10  *7 

1  : 

(Fpntijn  et  al.,  1964) 

(R51) 

Rates  given  at  300  °K  unless  otherwise  noted. 
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Reaction 

Reaction  Rate 
(cm3/sec) 

Reference 

N(2D)  —  N  +  hv  (5199A) 

T  =  103  sec 

(Wiese  et  al.  ,1966) 

(R52) 

O^S)  -  0(XD)  =  hv(5577A) 

t  =  0. 74  sec 

(Wiese  et  al.,  1966) 

(R53) 

0(XD)  -  0(3P)  +  hvgOAj 

t  =  110  sec 

(Wiese  et  al.,  1966) 

(R54) 

Competing  with  these  radiative 

transitions  are 

collisional  deactivation  reactions 

0(1S)  +  02  -  0*  +  0(3P) 

3  x  10-13 

(Zipf,  1969) 

(R55) 

0(1D)!+N2  -  N2  +0(3P) 

5  x  10  H 

(Zipf,  1969) 

(R56) 

0(1D)+02  -*02(a1Ag)  +0(3P) 

4  x  10'11 

(Zipf,  1969) 

(R57) 

and 

0(*S)  +0^0  +  0 

7.5  x  10~12 

(Felder  and  Young,  1972) 

(R58) 

For  large  NO  concentrations,  as  observed  during  auroral  activity 
(Zipf'et  al. ,  1970),  the  following  quenching  reactions  are  also  important: 


0(1S)  +  NO  -  (all  paths)  4  x  10~10  (Zipf,  1969)  (R59) 

and 

0(1D)  +  NO—  (all  paths)  1.5  x  10-10  (Zipf,  1969)  (R60) 


These  reactions,  however,  do  not  play  an  important  role  in  the  NASA 
electron  beam  experiment  where  the  ambient  atmosphere  was  quiescent  with  low 
NO  concentrations.  Three-body  reactions  are  very  slow  for  altitudes  above  100  km, 
and  are  not  included. 
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Diffusion  and  heating  can  play  an  important  part  in  the  chemical 
processes.  The  effect  of  diffusion  is  to  dilute  the  concentration  of  the  excited 
species  and  to  slow  down  the  reactions  between  two  excited  species.  The 
change  in  concentration  of  species  M  due  to  molecular  diffusion  in  Equation  (13) 
is  treated  as  an  axisymmetric  two-dimensional  problem.  The  diffusion  coefficient, 
D  ,  determined  by  Kenneshea  and  Zimmerman  (1970b),  is  shown  in  Figure  4. 

The  shaded  portion  of  the  curve  represents  the  eddy  diffusion  region  which  is 
generally  observed  below  106  km.  Above  this  altitude,  the  diffusion  becomes 
dominated  by  molecular  transport.  A  more  complicated  diffusion  coefficient 
affects  dissociative-recombination  processes,  Reactions  (R33)  through  (R37). 

Here,  the  diffusion  involves  two  charged  species.  Hence  coulomb  attraction 
and  the  effects  of  the  geomagnetic  field  must  also  be  taken  into  account. 


Electron  beam  heating  of 
the  irradiated  region  has  been  evalu¬ 
ated  by  accounting  for  the  energy 
deposition  when  one  ion  pair  is 
formed.  This  energy  is  distributed 
in  heating,  dissociation,  excitation, 
and  radiation.  A  detailed  evaluation 
of  heating  rates  requires  calculation 
of  the  electron  cooling  rate,  the 
chemical  reaction  rates,  and  the 
radiative  emission.  Since  this  is  a 
substantial  task,  the  heating  rate  is 
crudely  estimated  ior  the  atmospheric 
dose  levels  in  thn  NASA  experiments. 
The  energy  deposition  leading  to 
heating  the  atmosphere  is  estimated 
to  be  2/3  of  the  35  eV  required  for 


Figure  4.  Diffusion  Coefficient  at 

Different  Altitudes  (Taken 
From  Keneshea  and 
Zimmerman  (1970)) 
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ion  formation.  The  remaining  third  of  the  energy  is  assumed  to  be  either  used 
in  dissociation  or  radiated  away.  For  the  conditions  of  the  NASA  experiment, 
this  gives  a  temperature  rise  of  about  1°K.  Since  this  rise  is  small,  more 
detailed  calculations  have  not  been  considered. 
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SECTION  5 


APPLICATION  OF  THE  MODEL  TO  NASA 
"ARTIFICIAL AURORA"  EXPERIMENT 

In  this  section,  the  model  developed  in  the  previous  two  sections  is 
employed  to  analyze  th?  NASA  "artificial  aurora"  experiment.  We  shall 
determine  the  extent  of  the  electron  irradiated  region  and  compute  the  chemical 
and  radiative  processes  that  take  place  in  the  atmosphere . 

The  ionized  region  in  the  NASA  experiment  is  a  cylindrical  volume 
lying  along  the  geomagnetic  field.  Near  the  electron  gun,  the  beam  width  is 
nairow,  becoming  widest  in  the  region  of  maximum  energy  deposition.  The 
ion  pair  formation  rate  per  unit  volume  of  atmosphere,  calculated  from 
Equations  (2)  and  (10)  for  Eq  =  8.  7  kV  and  I  -  0.  49A  ,  is  shown  in 
Figure  5.  The  altitude  of  the  electron  beam  accelerator  was  taken  to  be 
250  km,  although  this  is  not  a  critical 
figure  as  very  little  energy  is  deposited 
above  200  km.  The  number  density, 
n(s)  ,  was  based  on  the  1962  Siandard 
Atmosphere. 

Since  each  electron  in  the  K 

UJ 

beam  slows  down  within  several  milli-  I 

seconds,  Lhe  atmosphere  is  effectively 

at  rest  during  the  energy  deposition  by 

an  electron.  The  total  ionization  per 

unit  volume,  which  depends  on  the 

irradiation  by  many  electrons,  however, 

does  depend  on  the  motion  of  the  elec-  FlgUre  5-  j°n  Production  Rate  With 

Respect  to  Altitude  Calcu- 

tron  accelerator  and  the  atmosphere,  lated  for  Conditions  of  the 

the  time  in  which  the  electron  gun  is  NASA  Experiment 

on  and  the  size  of  the  electron  cloud. 
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The  observed  size  of  the  irradiated  region  depends  also  upon  the  type 
of  observation.  That  is,  optical  measurements  of  prompt  emission  spectra 
would  indicate  a  small  cloud  size,  whose  lateral  spread  is  due  only  to  the 
electron  beam  multiple  scattering.  On  the  other  hand,  the  size  determined 
from  measuring  ti.e  emission  of  metastable  species  such  as  O^S)  ,  which  have 
long  lifetimes  compared  to  the  electron  beam  pulse,  would  appear  much  larger. 
Since  the  emission  from  0{  S)  persists  much  after  the  irradiation  time,  the 
observed  cloud  size  includes  not  only  the  lateral  spread  due  to  electron  scat¬ 
tering,  but  a  wake  of  emitting  air  left  by  the  moving  electron  beam. 

The  time  j  that  a  point  is  irradiated  for  a  continuously  operating 

beam  is 


T 


(17) 


where  is  the  lateral  spread  of  the  excited  cloud  (about  74  meters)  and  v^  is 
the  component  of  the  electron  gun  velocity  perpendicular  to  the  geomagnetic 
field.  The  latter  is  given  by  the  vector  equation 


x  (v  x  B)  _  v~-  (v* •  B)  B 


B 


B 


v-  (^)  cos 8  B 


(18) 


where 


cosfl  = 


IvjJD 

vB 


In  terms  Oi  .he  rocket  trajectory 


cos20  =  cos2  Y  -  [vH  cos{90  -  a  +  p  )  -  vv  sin(90  -  a  +  (3  )2| 


(19) 
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Here,  vH  and  vy  are  the  horizontal  and  vertical  components  of  the  vehicle 
velocity  calculated  from  the  apogee  and  range  of  the  rocket;  a  =  97°  is  the 
azimuth  of  the  rocket  trajectory  with  respect  to  true  north;  p  =  -8°  and 
y  -  69.5  are  the  declination  and  inclination  of  the  geomagnetic  field.  Note 
that  at  apogee,  where  vy  =  0,  Equation  (19)  reduces  to 

9  =  cos  1  sin (-15°)  cos  69.5°  =  95°  (20) 

The  velocities  across  field  lines,  v^  ,  for  the  four  visible  pulses 
(numoers  11  and  21  of  the  first  two  pulsi  lg  sequences)  are  presented  in  Table  3. 

TABLE  3 


VELOCITY  ACROSS  MAGNETIC  FIELD  LINES 


MEASURED  AT 

MEAN  BETWEEN 

CALCULATED 

FRANKLIN  CITY 

IGOR 

FRANKLIN  CITY 
AND  IGOR 

Pulse  1 

98 

68.6 

84 

110 

Pulse  2 

132 

134 

133 

131 

Pulse  3 

160 

187 

173 

154 

Pulse  4 

149 

224 

187 

178 

These  velocities  were  determined  from  observations  of  the  irradiated  region  from 
Igor  and  Franklin  City.  Notations  were  made  of  the  angular  velocity  of  the 
cloud,  the  distance  from  the  irradiated  region  to  the  observation  post,  and  the 
angle  made  between  the  line-of-sight  with  the  direction  of  motion  of  the  electron 
beam.  These  observed  velocities  are  presented  together  with  the  velocities 
using  Equation  (22)  with  the  horizontal  and  vertical  components  of  the  acceler¬ 
ator  velocity  der.  /ed  from  the  accelerator  trajectory.  A  sketch  of  the  motion 
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of  the  electron-beam-excited  cloud  is  shown  in  Figure  6.  The  irradiated  region 
is  approximately  clyindrical  in  shape  with  a  diameter  of  74  meters  and  a 
velocity  of  v  .  As  seen  in  Figure  6,  Region  A  is  illuminated  for  the  maximum 
time  that  the  atmosphere  can  be  irradiated  by  the  electron.  Neglecting  air 
motion,  this  time  corresponds  to  0.67,  0.56,  0.48,  and  0.41  sec  for  the  four 
brightest  pulses  respectively.  The  total  apparent  width  of  the  irradiated  air 
shown  by  the  distance  B  on  Figure  6  is 

S±(apparent)  =  S±+  v±rA  (21) 


where  rA  is  the  afterglow  duration  after  beam  shutdown,  0.38  sec.  If  winds 
are  again  neglected,  this  apparent  width  becomes  112,  123,  129,  and  141, 
respectively,  for  the  four  brightest  pulses.  Depending  on  wind  direction, 
atmospheric  winds  may  either  increase  or  decrease  this  irradiation  time. 

If  the  wind  is  blowing  in  the  same  direction  as  the  motion  of  the  electron  beam 
accelerator,  the  apparent  width  of  the  irradiated  cloud  becomes 


E  LECTROfJ  BEAM  ACCELERATOR  LOCATION 


S±(apparent)=  S±  +  (v±  -  VA)  ta  (22) 


Figure  6.  Sketch  of  the  Electron 

Irradiated  Region  as  the 
Electron  Accelerator 
Moves  Through  the 
Atmosphere 


where  vA  is  the  component  of  atmos¬ 
pheric  velocity  parallel  to  v  . 

The  atmospheric  emission 
produced  by  electron  beam  irradiation 
can  be  divided  into  two  oarts;  (1)  prompt 
(short  lifetime)  and  (2)  afterglow  (long 
lifetime)  radiation. 

The  prompt  radiation  in  the 
NASA  observations  was  from  the 
Ng(l-),  Ogfl-),  and  Ng(l+)  bands 
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as  produced  by  Reactions  (R3),  and  (R15),  respectively.  The  emission  from 
these  sources  is  shown  in  Figure  7  as  a  function  of  altitude.  Since  the  radiative 
lifetimes  of  these  sources  are  very  short  compared  to  the  electron  irradiation 
time,  their  radiation  remains  essentially  constant  over  the  electron  beam  pulse. 
The  upper  scale  in  Figure  7  gives  the  emission  integrated  through  the  irradi¬ 
ated  region  along  the  line-of-sight.  The  path  length  through  the  cloud  was  taken 
to  be  100  meters,  corresponding  to  line-of-sight  through  the  74  meter  wide 
cloud  at  108  km  altitude  when  seen  from  the  ground  at  Franklin  City.  The 
dashed  line  represents  the  combined  brightness  of  all  the  prompt  radiation 
sources,  and  the  dotted  line  represents  the  emission  from  these  radiators  when 
corrected  for  instrument  sensitivity.  Note  that  the  combined  brightness  of 
these  sources  exceeds  1  kilorayleigh  between  104  and  132  km.  This  is  in 
agreement  with  measured  results  presented  in  Table  1. 


Other  intense  prompt  radiation 
bands,  such  as  the  N2  second  positive 
and  the  Lyman-Birge-Hopfield  bands, 
are  also  excited  by  the  electrons.  For 
an  extensive  list  of  these  radiators,  see 
Chamberlain  (1961).  Only  the  radia¬ 
tion  which  falls  within  the  bandwidth  of 
the  NASA  detectors,  however,  are 
treated  in  this  section. 

The  radiation  from  O(^S) 
and  0(*D)  metastable  species  also 
contribute  to  the  observed  signal. 

The  intensities  of  these  radiators 
are  plotted  as  a  function  of  time  for 
106,  110,  and  120  km  in  Figures  8, 


EMISSION  OVER  THE  LINE  OF  SIGHT,  KR 


Figure  7.  Steady-State  Emission  Levels 
From  Prompt  Radiators  in 
the  3500  to  8000A  Region 
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Figure  8.  Calculated  Emission  Inten¬ 
sities  From  0(1S)  and  0(1D) 
During  Electron  Irradiation 
at  106  km  Altitude 


9,  and  10  respectively.  The  prompt 
3914A  signal  from  the  N^(1*)(00) 
transition  is  also  plotted  in  these 
figures  for  comparison. 

These  intensities  were 
calculated  from  the  mean  energy  depo¬ 
sition  late  at  each  altitude.  It  should 
be  noted  that  the  electron  irradiation 
is  not  uniform  throughout  the  cloud, 
but  instead,  is  more  concentrated 
near  the  beam  axis.  However,  even 
though  the  chemistry  changes  with 
the  variation  of  dose  rate,  the  total 
effect  is  somewhat  averaged  out  over 
the  width  of  the  cloud. 


Figure  9.  Calculated  Emission 
Intensities  From  0(!S) 
and  O(^D)  During  Electron 
Irradiation  at  110  km 
Altitude 


Figure  10.  Calculated  Emission 
Intensities  From  O(^S) 
and  O(lD)  During  Electron 
Irradiation  at  120  km 
Altitude 
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The  omission  from  O^S)  and  O^D)  is  seen  to  increase  with 
irradiation  time.  This  emission  may  continue  to  increase  even  after  the 
irradiation  due  to  the  continuation  of  dissociative- recombination  reactions 
which  are  a  source  of  these  states.  While  the  emission  from  directly  excited 
prompt  radiators  quickly  drops  to  zero  when  the  irradiation  ceases,  these 
metastable  states  continue  to  emit  ra^ition  for  several  seconds.  The  Of^Sj, 
which  accounts  for  a  major  part  of  this  radiation,  is  produced  in  Reactions 
(RIO),  (Rll),  (R17),  and  (R42),  and  is  depleted  by  radiation  (R53)  and  quenching 
(R58).  The  emission  in  the  afterglow  is  presented  in  Figures  11,  12,  and 
13,  for  altitudes  of  106,  110,  and  120  km,  respectively.  These  signals  were 
calculated  for  an  irradiation  time  of  0.67  sec,  corresponding  to  the  maxi¬ 
mum  irradiation  time  of  0.41  sec  were  also  carried  out  and  showed  the 
same  general  shape  with  only  slightly  lower  initial  intensities. 


time  tr  UK  imiouTm  »u.u. S((  *\*w  »t 


Figure  11.  Emission  Intensity  of  Figure  12.  Emission  Intensity  of 

Afterglow  Radiation  From  Afterglow  Radiation  From 

O^S)  and  O(lD)  After  0(*S)  and  0(*D)  After 

0.67  Sec  Irradiation  Time,  0.67  Sec  Irradiation  Time, 

106  km  Altitude  1 10  km  Altitude 
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for  the  direct  dissociative-excitation 
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ital  with  the  same  intensity  aa  that 
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follcmtng  the  1I17A  ultraviolet  emission  from  dissociative -excitation  Into  the 
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SECTION  6 


01 SCUSS ION  AND  CONCLUSIONS 

Even  though  the  data  obtained  from  the  NASA  Artificial  Aurora 
experiment  were  limited,  the)'  allowed  us  to  derive  several  Important  results 
concerning  the  effect  of  electron  bombardment  of  the  upper  atmosphere.  The 
sue,  shape,  and  altitude  dependence  of  the  observed  signal  gave  some  Insight 
to  how  electrons  deposit  their  energy  Into  the  atmosphere  nnd  the  upper  atmos- 
phrrle  chemistry  involved  in  a  disturbance  by  energetic  electrons. 

The  electron  slowdown  In  the  upper  atmosphere  has  been  modeled 
by  a  simple  attenuation  process  by  fitting  Equation  (1)  to  the  altitude  dependent 
atmo»|»hcrlc  number  density.  The  site  and  shape  of  the  cloud  has  also  been 
modeled.  The  electron  interaction  with  the  atmosphere  was  modeled  by  fit¬ 
ting  the  energy  distribution  as  measured  In  the  laboratory  'o  a  varying  density. 
The  radial  spread  of  the  Irradiated  region  was  calculated  by  replacing  the 
electron  mean -free -pslh  by  the  electron  Larmor  radius  The  resulting  electron 
interaction  using  this  lechnique  was  consisiant  with  the  results  of  Berger 
et  al.  (1970). 

In  the  NASA  experiment,  the  site  and  shape  ol  the  cloud  were  deter¬ 
mined  from  measurements  o.‘  visible  light  emission.  In  order  to  Interpret  the 
data,  both  the  Interaction  of  electrons  with  the  atmosphere  and  the  resulting 
chemistry  and  radiation  have  been  calculated.  This  visible  emission  consisted 
of  essentially  two  classes  of  signals:  prompt  (short  lifetime)  radiation  which 
disappears  in  less  than  a  millisecond  after  beam  turn-off,  and  long-llfetlmc 
radiation  from  metastablc  states  whose  radiation  persisted  for  several  tenths 
of  a  second  after  beam  turn-off.  The  width  of  the  observed  column  has  been 
determined  by  combining  the  radiation  temporal  distribution,  the  Instrument 
response,  and  the  motion  of  the  electron  beam  accelerator.  While  the  data 
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taken  during  the  NASA  experiment  shows  too  much  scatter  to  infer  an  exact 
cloud  width,  the  calculated  cloud  sizes  agree  with  the  experimental  observations 
within  experimental  error. 

The  motion  of  the  electron  beam  accelerator  determined  the  air 
irradiation  time.  Atmospheric  winds,  which  have  not  been  considered,  would 
move  neutral  species  in  and  out  of  the  irradiated  region  and  complicate  the 
analysis.  The  effect  of  winds  on  neutral  species  is  to  increase  or  decrease 
their  interaction  time  with  the  electron  beam.  The  electrons,  on  the  other 
hand,  are  trapped  by  the  magnetic  field,  and  are  affected  by  the  atmospheric 
motion  only  through  their  attraction  to  the  ions.  The  ions  are  dragged  by 
collision  with  neutral  species,  and  by  coulomb  attraction  to  the  electrons. 

Thus,  including  atmospheric  motion  in  the  analysis  would  change  , 

the  afterglow  signal  due  to  metastables,  while  the  prompt  radiation  would 
remain  essentially  unchanged. 

The  chemistry  involved  in  the  electron  beam  excitation  yields  inter¬ 
esting  results  even  from  the  limited  data  available.  Detailed  study  of  the  after¬ 
glow  signal  from  the  NASA  Artificial  Aurora  determined  the  amount  of  0(*S) 
produced.  While  the  prompt  radiation  signal  consists  ot  the  sum  of 

N„(l+),  and  O-fty)  bands,  the  main  contributor  to  the  afterglow  is  C 

‘  ‘  2  1 
Other  metastable  radiators  in  the  visible  include  N(  D)  and  0(  D)  which  have 

much  longer  radiative  lifetimes  and  thus  contribute  insignificantly  to  the  total 
emitted  signal.  Analyzing  the  dv.ta  available,  as  discussed  in  the  previous 
section,  leads  to  the  conclusion  that  the  O(^S)  production  can  be  attributed  pri¬ 
marily  to  the  dissociative-recombination  of  ot  . 

As  seen  from  Figures  8,  9,  and  10,  the  total  0(1S),  5577A  signal  nay 
build  up  to  over  one  third  of  the  Ng(l-)  3914A  signal  within  one  second  of  irradia¬ 
tion,  in  agreement  with  what  has  been  observed  in  auroras.  The  additional 
contributions  of  the  (v'  =  0,  v"  =  1)  and  higher  vibrational  bands  of  N„(l+)  bands, 

as  well  as  the  Q„(l-)  band,  however,  make  the  0(  S)  signal  much  lower  than  the 
4  1 

prompt  radiation  signal.  This  0(  S)  signal,  when  corrected  for  instrument 
sensitivity,  is  consistent  with  the  NASA  observations. 
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The  model  developed  in  this  report  for  calculating  the  visible 
emission  from  electron  beam  excitation  shows  good  agreement  with  observed 
results.  While  these  calculations  were  used  to  analyze  the  NASA  experiment 
performed  in  1969,  the  same  techniques  should  be  valid  for  predicting  the 
results  of  similar  types  of  atmospheric  disturbances  including  future  electron 
beam  experiments. 
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